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ABSTRACT 

^! We investigate the effect of planetary rotation on the transit spectrum of an 

2 '■ extrasolar giant planet. During ingress and egress, absorption features arising 

c/3 , from the planet's atmosphere are Doppler shifted by of order the planet's rota- 

tional velocity (~ 1 — 2 km s~^) relative to where they would be if the planet were 
^ ! not rotating. We focus in particular on the case of HD209458b, which ought to be 

^ ! at least as good a target as any other known transiting planet. For HD209458b, 

O ■ this shift should give rise to a small net centroid shift of ~ 60 cm s~^ on the stel- 



O 
p 

in 



lar absorption lines. Using a detailed model of the transmission spectrum due 
to a rotating star transited by a rotating planet with an isothermal atmosphere, 
O I we simulate the effect of the planet's rotation on the shape of the spectral lines, 

Q I and in particular on the magnitude of their width and centroid shift. We then 

use this simulation to determine the expected signal-to-noise ratio for distin- 
I guishing a rotating from a non-rotating planet, and assess how this S/N scales 

^ I with various parameters of HD209458b. We find that with a 6 m telescope, an 

equatorial rotational velocity of ~ 2 km s~^ could be detected with a S/N ~ 5 
by accumulating the signal over many transits over the course of several years. 
With a 30 m telescope, the time required to make such a detection reduces to 
less than 2 months. 

Subject headings: astrobiology - planetary systems - radiative transfer - stars: 
atmospheres - stars:individual (HD209458) - astrochemistry 



1. Introduction 

Since the early 1990s, more than 200 planets have been discovered orbiting other stars. 
Roughly four fifths of the planets that have been discovered are at least half as mas- 
sive as Jupiter, and about a quarter of the known planets orbit extremely close to their 
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parent star ( ^ 0.1 AU). The Jupiter-mass planets that are in close orbits are highly 
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Observations of the transiting planets have confirmed their similarity to Jupiter by revealing 
that these are gas giant planets, with radii comparable to, or somewhat larger than Jupiter's 
Jcharbonneau et al.l[200ol : Iflenrv et al.lboool : lGaudill2005h . 



The effects of the tidal torques experienc ed by an orbiting body h ave been studied for 
a long time - for an early seminal analyis, see iGoldreich &: Peald ( 119661 ). Such torques tend 
to synchronize a satellite's rotation rate to its orbital rate, and if the torque is sufficient 
this synchronization is achieved and the orbiter is said to be "tidally locked," as the Earth's 
Moon is. The Hot Jupiter-class extrasolar planets are thought to orbit sufficiently close to 
their stars that their tidal locking timescales are much shorter than the ages of the planets. 
The planets, then, are expected to be tidally loc ked to the stars, with o ne hemisphere in 
permanent day and the other in permanent night ( Harrington et allbood ). 



A tidally locked Hot Jupiter will have a permanent sharp contrast in temperature 
between the substellar point a nd the night side, which must have a profound influence 
on the atmospheric dynamics. IShowman fc GuillotI (|2002| ) make simple predictions of the 
day/night temperature difference (~ 500 K) and the speed of winds (up to ~ 2 km s~^), 
and their detail ed, three-d i mens ional simulations agree with their estimates. Shallow-water 
simulations by ICho et al.l (120031 ) predict longitudinally averaged zonal wind speeds of up 
to 400 r n s~^, with local winds app roaching 2.7 km s~^ (under some assumtions). Simula- 
tions by ICooper fc ShowmanI (120051 ) predict a super-rotational jet (i.e., blowing eastward, 
where north is deflned by the right-hand rule) that blows the hottest part of the planet 
downstream by about 60° from the substellar point. Their simulations predict supersonic 
winds exceeding 9 km s~^ at high latitudes, high in the atmosphere (where the optical depth 
is low) and winds exceeding 4 km s~^ at pressures near the photosphere. A Spitzer Space 
Telescope phase curve for v Andromedae b rules out a phase-sh ift as large as 60° between 



200' 



^, but a Spit zer phase curve 



the substellar point and the hottest spot (iHarrington et al. 

for HD189733b favors a ~ 30° shift for that planet (iKnutson et al.l l2007al ). so it remains 



unclear to what extent available data indicate very strong photospheric winds. 



Transmission spectroscopy is a way to probe the atmospheres of these planets. ICharbonneau et al 



( I2OO2I ) were the flrst to detect an absorption feature in what is probably the atmosphere of 
HD209458b, when they found that the effective radius of the planet increases slightly at the 
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waveleng th of a strong sodium absorptio n doublet (the sodium D lines) at ~ 590 nm. In 
addition, IVidal-Madjar et al.l (l2003l . 1200J) have reported a number of absorption features in 
HD209458's transit spectra that are due to various species (including hydrogen Lyman al- 
pha, neutral carbon, oxygen, and sulfur, and some ionization states of carbon, nitrogen, and 
silicon) in a hot exosphere that is probably only loosely bound to the planet. Intr iguingly, 
throug h analyzing the red and near-IR portion of HD209458b's transit spectrum iBarman 
( 120071 ) found a lOa detection of atmospheric water vapor. Several measurements of the 
planet's emission spectrum, however, have fou nd results that seem to be inconsistent with 
high water abund ance high in the atmosphere (IGrillmair et al.ll2007l : iRichardson et al.ll2007l : 
Swain et al.ll2007h . 



Initial work by ISeager &: Sasselovi ( l2000l ) and a comprehensive study by iBrownl ( 12001 



hereafter BOl) have described various other considerations that should affect the details 
of transit spectra, including the orbital motion of a planet (a few tens of kilometers per 
second in the radial direction), the rotation of the planet (a few kilometers per second at 
the equator, according to the hypothesis that the planet is tidally locked), and winds on the 
planet's surface (in BOl's analysis, up to ~ 1 km s~^). These physical effects should tend 
to broaden or impose Doppler shifts on absorption features due to the planet's atmosphere. 
BOl constructed an impressively detailed model of radiative transfer through a Hot Jupiter's 
atmosphere, assuming various models of zonal windfiow superimposed on an equatorial bulk 
rotation speed of Veq = 2 km s~^, which is approximately the value for HD209458b under 
the assumption that it is tidally locked in its 3.5 day orbit. He finds the height of the cloud 
deck to be the most important parameter that affects the transmission of light through the 
planet's atmosphere. 

The original discovery of the rough ly Jupiter-mass planet in a close, ~ 4 day orbit 



around 51 Pegasi (IMayor &: Queloall995l ) prompted i nterest in the dynam ics and structure 



that must govern a highly insolated gas giant planet (iGuillot et al.l 119961 ). Observations of 
the transiting Hot Jupiters heightened this interest when they revealed a puzzling feature of 
these planets: at least several of them are a bit puffier than Jupiter, with diameters ranging 
from slightly larger than Jupiter's to as much as ~ 80% larger. It is not clear what allows 
some planets to maintain such large radii. It has been suggested that if a Jovian planet 
migrates very quickly, from its presumed formation location at least several AU from its 
star, to its eventual several day orbit, then it might reach its final home before it has cooled 
enough to shrink to Jupiter's radius. Accordingly, some authors have investigated the mi- 
gration processes t hat lead gas giant planets to such close orbits as have been found (e.g. 



Trilling et al.ll2002l ). Others have investigated various ways in which a gas giant could either 



sustain its infiated size ( 


Guillot & Showman 


2002; 


Burrows e 


t al. 


2003; 


Laughlin et al. 


2005; 


Bodenheimer et al. 


2003: 


GuillotI 


2005; 


Burrows et al. 


2007; 


Chabrier & Baraffe 


2007 


). Al- 



though various physical mechanisms have been suggested as the apparently missing energy 
source that allows the unexpectedly large radii sometimes seen, the lesson of these investiga- 
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tions in toto is that it is not easy to explain the inflated sizes, either in terms of the greater 
stellar flux that these planets experience by virtue of being so cl ose to their stars , or in 
terms of their evolutionary migratory histories. A recent paper by IWinn fc Holmaru (120051 ) 
propose that, contrary to the commonly accepted paradigm. Hot Jup iters might be trapped 
in a Cassini state with large obliquity, in which the spin-axis precesses in resonance with the 
orbit, but lies nearly in the orbital plane. Such a state might be stable against perturbation, 
and yet able to generate sufficient internal energy to increase a gas giant planet's radiu s 
to the observed values. In light of an even more recent analysis by iLevrard et al.l (120071 ). 
however, it appears that the probability of capture into a Cassini st ate 2 resonance i s quite 
small for a planet with semi-major axis a < 0.1 AU. Furthermore, iFabrycky et al.l (120071 ) 
argue that even if a planet is captured into Cassini state 2, it is likely to remain there for a 
time that is short relative to the age of the system. 

High-resolution transit spectra that have high signal-to-noise ratios will allow us to 
distinguish between various models of orbit, rotation, and weather, as discussed by BOl. 
Because the orbit is known to high accuracy, and the predictions of the effects of weather 
(or climate) are highly uncertain, as described above, we will focus in this paper on the 
much more easily predicted effect of a planet's rotation on a transit-spectrurn . If we neglect 
winds, then the large-obliquity Cassini state described by IWinn fc HolmanI (120051 ) should 
have a spectral signature that is very similar to that of a non-rotating model. In contrast, 
the rotation of a tidally locked planet should impose a Doppler distortion on spectral lines 
arising from the planet's atmosphere that is roughly an overall redshift during ingress, as 
the planet is just entering the stellar disk, and a similar distortion that is roughly an over- 
all blueshift during egress, as the planet is just exiting the disk. During mid-transit, the 
spectral distortion is more similar to rotational broadening. In the present investigation, 
we address whether there is any hope that these spectral distortions from tidally-locked 
rotation can be obs e rved. In our study, we focus only on the sodium doublet detected by 
Charbonneau et al.l (120021 ). As we will show below, the sensitivity of a measurement of ro- 
tation scal es with the square roo t of t he number of li nes under consideration. Model spectra 
from, e.g.. ISudarsky et al.l (120031 ) and iBarmanI (120071 ) predict a strong potassium doublet at 
~ 770 nm, strong water absorption features in the near-infrared, and a handful of near-UV 
lines. If some of these are confirmed in the atmosphere of a transiting planet, they will 
provide a modest increase in S/N. Since the sodium lines are expected to be the strongest, 
however, it seems unlikely that observing multiple lines will yield a boost in S/N by more 
than a factor of a few. 



We emphasize that it may not be at all justified to neglect winds. It is quite likely that 
there are super-rotational winds on Hot Jupiters, which are probably necessary to heat the 
"night" side. As indicated above, some models predict, and the observed phase curve for 
HD189733b suggests, that at the photosphere these winds might be significantly (100% or 
more) greater than the equatorial rotation rate, and therefore might contribute importantly 
to the Doppler distortion induced by the motion of the planet's atmosphere. Nevertheless, 
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in order to isolate the contribution of rotation, we do neglect winds in this study. The 
Doppler distortions that we predict can therefore probably be taken as a lower bound on the 
distortions that would be observed for a tidally-locked transiting Hot Jupiter. 

We find that the spectral shifts induced by rotation will be difficult to detect with 
current technology, but perhaps not insurmountably so, at least with technology that might 
be available in the not-to-distant future. The measurements we will describe are limited by 
a paucity of photons. As such, their signal-to-noise ratio will be enhanced by a bright star 
and a puffy planet (i.e., a planet with a large scale-height). HD209458 is at least a magnitude 
brighter than any other star with a known transiting planet except HD 189733, and its planet 
is larger than HD189733b; so HD209458b should be a better target than any other known 
transiting planet except possibly HD189733b. In this paper, we model the HD209458b 
system because it is the best-studied system, and it is unlikely that any currently-known 
planets would be significantly better targets. In a single transit, observations of HD209458 
with a 6 m telescope that has a high-resolution (> 50, 000) optical spectrograph with good 
throughput (~ 18%) could only show the infiuence of tidally locked rotation at the ~ 0.2a 
level. With ultrahigh-resolution ( ^ 700, 000) and good throughput (~ 4%) this effect would 
still only show up at the ~ O.Gcr level. In less than a year, the signal of rotation could be 
present at five times the noise (5'/A^ = 5). Of course, a telescope with larger collecting area, 
higher spectral resolution, or better throughput would cause the signal to be apparent at 
that significance level in less time. 

Other stud ies have approached the problem of determining t he rotation rate from a 



different angle. ISeager &: Hull (120021 ) and iBarnes &: Fortneyl (120031 ) suggest that an oblate 



spheroid will have a different transit light curve from a perfect sphere, and so measuring the 
oblateness from transit photometry will provide a handle on the rotation rate. The oblateness 
is somewhat degenerate with several other parameters that are not perfectly known, however, 
so they conclude that it would be difficult to actually determine that rotation rate in this 
manner. The method we describe here could eventually prove to be an important complement 
to other observations to constrain the rotation rate. 

In the remainder of this paper, we address this idea in detail. One complication that 
we discuss below is that the technique of this paper is not immune from several near- 
degeneracies among the many attributes of transiting extrasolar planets that infiuence light 
curves or spectra. Although it is likely that current or near-future instruments will be 
sensitive enough that the spectral distortion imposed by HD209458b's rotation (if it is tidally 
locked) is visible, it might still be very challenging to discern the fingerprint of rotation from 
other attributes that affect the spectra at a similar level. In this paper, we tackle the 
forward problem of calculating the amount of distortion that is caused by rotation. The 
inverse problem - determining from observations whether a planet is tidally locked - is more 
difficult and should be the topic of a future study. 

The structure of the rest of this paper is as follows: In § [2l we describe qualitatively what 
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happens to the starhght received on Earth when a planet transits its star; we give a rough 
order of magnitude estimate of the the magnitude and detectabihty of the spectral distortions 
caused by tidally locked rotation; and we briefly describe some technological progress and 
remaining challenges relevant to our task of deducing bulk motions in a planet's atmosphere 
from transit spectra. In § [3], we describe our computational model of a transit spectrum. In 
§ m we describe the results of our model according to various assumed input parameters. In 
§13 we discuss the scaling of S/N on various model parameters and we address the prospects 
of actually observationally determining whether a transitting planet is tidally locked. In § |6l 
we conclude by describing various ways to boost our predicted S/N to a more optimistic 
value. 



2. Overview of the Problem 

The practical feasibility of the investigation we undertake depends on a few factors: 
understanding the various detailed processes that affect the starlight that reaches Earth 
when a planet transits its star; the magnitude of the distortion that tidally locked rotation 
induces; and the technology available to measure such distortions. In this section, we give 
an overview of these three factors - in particular, in § \'2.2\ we give a simple estimate of the 
results that we will later (in § Hj) calculate in detail. 



2.1. Relevant Processes 

A planet transiting in front of its star affects the starlight that ultimately reaches Earth 
in many ways. The motion of the planet's atmosphere (rotation and winds) is a small 
perturbation on top of several more dominant effects. We therefore summarize below the 
physical processes that are at least as significant as the effect of tidally locked rotation. 
Figure [1] schematically represents this situation, and captures nearly all of the processes 
described below: a rotating planet (of exaggerated relative size) transits in front of a rotating 
star. The figure depicts a snapshot partway through ingress, when half of the planet is in front 
of the stellar dislo- The white circle indicates a hypothetical sharp demarcation between 
the opaque part of the planet (in black) and the optically thin part, labeled "Atmosphere" 
(described further below). 

1. Geometric Occultation: 

The largest effect is an overall dimming by a factor of roughly the ratio of the area of 
the planet to that of the star: {Rp/R^Y. Since stars are not perfectly uniform disks. 



^The planet is above the star's equator to represent a slight inclination in its orbit. 
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Fig. 1. — Rotating planet beginning to transit in front of rotating star. The vertical black 
arrows represent the rotation axes of the planet and the star, and the curved arrows indicate 
the direction of rotation for each. The X's on the right-sides of both the planet and the 
star indicate regions that are receding from the observer and are therefore redshifted; the 
unmarked left sides of the planet and the star are moving toward the observer and are 
therefore blueshifted. The white circle surrounding the opaque black part of the planet 
denotes the cloud deck, or the boundary between the partially transparent and the fully 
opaque portions of the planet's disk. The planet is orbiting in the same sense as both it 
and the star are rotating. The planet is shown above the star's midplane to represent the 
inclination of the orbit relative to the line-of-sight. 

but instead tend to darken toward the limb at most visible wavelengths, the fractional 
dimming due to being in the planet's shadow tends to be slightly less than the ratio 
of the areas when the planet is near the edge of the stellar disk and slightly more than 
this ratio when the planet is near the center. 

2. Stellar Wobble: 

The primary spectral effect of the planet orbiting the star is the radial velocity wobble 
induced by the planet's gravity. This periodic spectral shift is of course in effect during 
the transit, when, for a close-in planet like HD209458b, it has an influence on the order 
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of ~ ±10 m s ^. This effect is a redshift as the planet begins to transit across the disk 
(during ingress) and a blueshift during egress. 

3. Rossiter-McLaughlin Effect 

A more subtle effect arises because, during the transit, the planet moves across - and 
therefore blocks - parts of the star that have different recessional velocities. If (as is 
expected) the planet's orbit is aligned with the star's spin, then during ingress the 
planet is blocking a part of the star that is slightly blueshifted, and during egress it is 
blocking a part of the star that is slightly redshifted. Figure [1] illustrates the planet 
blocking some of the bluest parts of the star during ingress. 

The parts of the star that are occluded during ingress/egress have spectra that are 
blue/redshifted by a velocity that is approximately the equatorial rotational speed of 
the star, or about ~ 1-2 km s~^ for a Sun-like star. As the figure indicates, during 
ingress/egress, the integrated spectrum of the remaining (unblocked) parts of the star is 
on average slightly redder/bluer than it would be if the planet were entirely transparent. 
Therefore, during ingress, the centroids of stellar lines are shifted slightly to the red, 
and during egress the centroids are correspondingly shifted to the blue. 
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already caused by the radial velocity induced by the planet's gravity, described in 
(2.) above. The R ME has beeri descr ib ed in depth mo re rec ently in the context o f 
extrasolar planets bv lohta et aP J2005h . Icim^nej J2006h . and ICaudi fc WinnI J2007h . 



These centroid-shifts are expected to be comparable in magnitude to the radial velocity 
wobble from the planet's gravity, and can be roughly estimated as 



M 



1 km s"^ X {Rp/R^Y ~ 10 m s" 



In fact, the amount of the shift can be predicted precisely for a given orientation of the 
planet's orbit, and so measuring the shift is tantamount to measuring the alignment 
between the star's spin and the planet's orbit. Three years ago. I Winn et al.l (120051 ) first 
found that the spin of IID209458 and the orbital plane of its planet are nearly a l igned . 
The degree of alignment has been measured for two other systems - IWinn et al. J2006h 



found that the spin of HD189733 and its planet are also nearly aligned, and lNarita et al. 



( 120071 ) measured a mis-alignment between these two vectors by ~ (30 ± 20)° in the 
TrES-1 system. 

4. Planet's Atmospheric Opacity: 

Furthermore, a gas-giant planet's opacity surely does not have a perfectly sharp dis- 
continuity at an outer boundary, as a billiard ball does. Instead, it has an extended 
atmosphere in which the opacity must vary more or less smoothly. There may be 
a cloud layer, below which the planet is entirely opaque to tangential rays and above 
which the opacity varies smoothly. Most critical to our investigation, at a given radius, 
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the planet's opacity to tangential lines of sight must vary with wavelength, depending 
on the contents of its atmosphere. At wavelengths of strong atomic or molecular tran- 
sitions, the planet's atmosphere will be more opaque than at other wavelengths. As a 
result, the effective radius of the planet, or the radius at which the optical depth along 
a tangential ray is of order unity, is greater at some wavelengths than at others. These 
effects have been described in detail by BOl. 

5. Planet's Orbital Motion: 

The motion of the planet's atmosphere must influence the transit spectrum in several 
delicate ways. As BOl points out, there are three main mechanisms by which the 
motion of a planet's atmosphere relative to its star can affect the spectrum: the planet's 
orbital velocity along the line-of-sight, the planet's (possibly tidally locked) rotation, 
and winds in its atmosphere. The largest effect of these three is the orbital velocity, 
which imposes a bulk blue/redshift during ingress/egress of ~ 15 km s^^ to spectral 
lines arising from the planet's atmosphere. These shifts are of opposite sign to the 
radial velocity wobble and to the shifts from the RME, and therefore tend to lessen 
the apparent RME slightly. 

6. Planet's Atmospheric Motion: 

The most djTiamically interesting (and subtlest) effects are those caused by the plane- 
tary rotational velocity and atmospheric winds. Since a tidally locked planet rotates in 
the same sense as it orbits, the rotational velocity of its outside edge has the same sign 
as its orbital velocity, and the rotational velocity of its inside edge has the opposite 
sign. As a result, during the beginning of ingress and the end of egress, when only the 
inside edge of the planet is over the star, tidally locked rotation will impose a spec- 
tral distortion that is in the opposite sense of that caused by the bulk orbital velocity 
described in (5.) above, and that is in the same sense as the RME: the distortions 
are roughly equivalent to a relative redshift during ingress (graphically represented in 
Figured]) and a relative blueshift during egress. During mid-transit, with some parts 
of the rotating planet's atmosphere moving toward and other parts away from the star 
relative to an otherwise identical but non-rotating planet, the overall influence of the 
planet's rotation is approximately equivalent to rotational broadening. 

Winds complicate the picture even further. It is likely that winds tend to rush from 
the substellar hot spot to the colder night side of the planet. With the substellar point 
on the opposite side of the planet from Earth during a transit, this corresponds to 
winds rushing toward us at several hundred to several thousand meters per second. 
This would tend to blueshift the spectrum throughout the transit. Zonal wind bands, 
somewhat similar to those on Jupiter but with much higher speeds, or other more 
detailed winds, can have an even more intricate effect. 

7. Additional Effects: 

If a transiting planet were to have nonzero orbital eccentricity, or rings, these could 
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complicate a measurement of rotation rate. Nonzero eccentricity would break the sym- 
metry between ingress and egress. Still, if the orbit were well-known, this could be mod- 
eled and taken into account. It seems unlikely that a Hot Jupiter could maintain rings: 
Icy rings would sublimate, and, if not continuously repl enished, dusty / r ocky rings 
would quickly succumb to the Poynting-Robertson effect fjPoyntingl 1 19031 : [Robertson 



19371 ). But if, somehow, a ring were to find a way to persevere around a transiting 



Hot Jupiter, it could confound - perhaps hopelessly - a measurement of ro tation. The 



conseq uences of rings for the Rossiter-McLaughlin effect is addressed in lOhta et al. 



(120061 ). Saturn's rings are nearly four-times the area of the planet, so for a planet 
(with equatorial rings that are as relatively large as Saturn's) whose orbit is tilted an 
angle 0.1 (in radians) from edge-on, the rings would be ~ 40% the area of the planet, 
which would increase the RME by ~ 40%. Uncertainty about the presence and size of 
a ring introduces an uncertainty in the size of the RME effect that is probably larger 
than the size of the rotation effect. Furthermore, a ring would occlude a (small) part 
of the planet's atmosphere, which would (slightly) reduce the strength of the rotation 
signal. 



Other interesting phenomena that prima rily affect a transit lig ht curve, rather tha n 
the spectrum, include star-spots Jsilva 2003 ). atmosph eric lensing (Hui fc Seager 2002), 
and finite-speed-of-light effects ( Loebl 20051 ). Although Winn fc Holman ( 2005 ) describe 
a possible configuration (Cassini state 2) that would produce a spectral signature that is 
nearly identical to what would be expected from a non-rotating planet, the likelihood that 
any Hot Jupiters are in such a configuration might be low, and it seems quite likely that 
some transiting planets are not in this state. Nonetheless, the motion of a transiting planet's 
atmosphere - rotational, wind, or other - is clearly interesting, and the basic technique that 
we describe below is applicable to any model of atmospheric motion. 



2.2. Preview of Results 

A rough estimate of the velocity-shift that is imposed during ingress to the centroids 
of the stellar Na D-lines by the planet's tidally locked rotation (on top of the RME and the 
shift from the planet's orbital velocity, both of which would be present even if the planet 
were not rotating) is the following: 




~ 0.64 X 1% X 15% X 2000 m s" 



= 1.9 ms~^ (1) 

In this equation, is a planet-centered azimuthal angle, Rp and R^, are the planet's and star's 
radius, respectively, Hatm is the height of the planet's atmosphere, and Veq is the equatorial 
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rotation speed. The rotation speed at angle is v^q cos[0]. We take the average of cos[0] from 
—tt/2 to n/2 to get the average planetary rotation speed. We have used Ilatm = 7500 km, 
or 15 times the presumed scale height of 500 km, because the sodium lines are so heavily 
saturated that at the assumed abundance and cl oud deck height in our mod el the lin e cores 
do not become optically thin until that height. iBurrows et al.l (120041 ) and iFortneyl (120051 ) 
describe how the optical depth along tangential rays is greater than the optical depth on 
normal rays. The product 



atm 



2 -R* 



STT-RpIIatm 



Rp 



n 



atm 



Rn 



is the ratio of the area of the portion of the planet's atmosphere that is in front of the star 
halfway through ingress to the total area of the disk of the star. Based on this estimate, 
we expect a maximum velocity shift of 6v ~ 190 cm s~^. If we take into account that 
HD209458b's orbit is actually slightly inclined relative to the line of sight, the cosine average 
decreases to ~ 0.45, and the total estimate decreases to ~ 140 cm s~^ This estimate is in 
reasonably good agreement with the centroid-shifts predicted by the full model calculation 
below (~ 60 cm s~^); the difference between the estimates is most likely due to the difference 
between the shapes of the stellar and planetary lines. 

We now estimate the signal-to-noise ratio for the detectability of this effect in an ob- 
servation of duration At, with a telescope that has diameter D and throughput efficiency 
7]. The signal is the distortion of the spectrum relative to a non-rotating planet, and for 
now we will assume that the noise is dominated by photon noise. If a spectrum F[X\ with a 
symmetric absorption feature of depth AF centered at Aq is redshifted by an amount AA to 
F[X] = F[X — AA], what is the integrated absolute difference \F — F\ over some wavelength 
range 2L centered on Aq? If the absorption feature is wide compared with AA, then, by 
symmetry. 



S 



[' \F[X]~F[X-AX]\dX 
2 [ °^ {F[X] - F[X - AX]) dX; 



(2) 



and if AA is small then 



/•Ao+L 

S ^ 2AX F'[X]dX 
^ 2(AA)(AF). 



(3) 



We may now estimate the S/N of our effect (for a single absorption line) using the 
lesson of equation ([3]), provided we know the absolute normalization of the stellar spectrum 
(the number of photons per unit wavelength). A spherical blackbody of radius Rj, and 
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temperature T^,, at distance d from the telescope, has a photon number flux at wavelength 
A of 



dX ' ' ' \hc/X J \ d^ 



2 



Bx[\T] -— xr^n{D/2) 



> (4) 



2A4(exp[(/ic)/(AA;T,)] - 1) ' \ d 

where B\ is the Planck function. Since the fractional decrease in the spectrum at the 
line-center is approximately 5atm; we may express the parameter AF from equation ([3]) as 
AF 6atm{dN^//dX). Similarly, since the root-mean-square velocity shift during ingress 
is (f^)^^^ ~ (1/2) X (2000 m s^^) = 1000 m s^^jfl we may express the parameter AA as 
AA ~ ((^^) /c) X Aq. The distortion (the signal) from a single line can therefore be 
estimated as 

S = SN.r^ 2(AA) (^S^tm^^ At 

X V (At). (5) 



A4 {exp[{hc)/{Xkn)] - 1) ' \ d 



The shot-noise is the square root of the number of photons in a wavelength range 2L roughly 
equal to the FWHM of the line, or about 7 km s~^ for a heavily saturated line such as the 
Na D lines under consideration: 



N ~ ^/^(2L)(At) 



AMexp[(/ic)/(AA;r,)] - 1) ^ / 

We estimate the total signal-to-noise ratio arising from a single absorption line, during an 
ingress integration of duration At, to be roughly 

~ (6.6 X 10"^) (3.3 X 10"^) (2.1 x lO") (5.0 x 10"^) ^ 

~ 2.3^/^. (7) 

The above calculation uses parameters for HD209458 and its planet, a sodium D line, and 
a 6 m telescope: A = 600 nm; At = 1000 s; R, = 7.3 x 10^° cm; T, = 6100 K, d = 47 pc; 
and D = 600 cm. For two identical absorption lines, we gain a factor of 2^/^ in S/N, and for 



•^We write (1/2) x (2000 m s ^) because the mean value of cos^ from —tt/2 to tt/2 is 1/2 
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egress we gain another factor of 2^/^, giving a total one-transit S/N of roughly 4.67]^/^, not 
counting the additional signal available during mid-transit (see further discussion below). 
This S/N ratio is in principle independent of the spectral resolution of the spectrograph, for 
sufficiently high spectral resolution. For low spectral resolution, however, the S/N could be 
lower than this estimate (below, we conclude that the S/N loses its dependence on resolving 
power for spectral resolution ^ 500,000). 

There were several optimistic assumptions that went into this estimate. Still, this rough 
estimate of the degree to which a planet's rotation influences its transit spectrum indicates 
that the more in-depth study that we perform below is warranted. 



2.3. Available Technology 



Detecting the centroid-shifts caused by tidally locked rotation ( ^ 1 m s~^) will require 
very precise measurements of stellar transit spectra. Obtaining such high precision spectra 
will be quite challenging, for a number of re asons, several of which were described in the 
groundbreaking paper by iButler et al.l (119961 ) that analyzes the limits of Doppler precision. 
Of particular concern, stellar pulsations and turbulent motions in ste llar photospheres can 
cause small regions of the stellar disk to move at up to 300 m s~^ (IDravinsi Il985l : lUlrich 
199ll ). These motions tend to average out to produce stellar spectra that are largely stable; 
but it is likely that at least some giant convection cells are not small relative to the size of 
a planet, and these could introduce a contaminatin g sour ce of noise when they are located 
behind the planet or its atmosphere. iButler et al.l (119961 ) reviewed what was then known 
about the variabihty of stellar line-proflles; the upshot is that line-widths may vary by up 
to several meters per second over several years, but it is not clear to what extent spurious 
apparent velocity shifts may be induced by convection, and such stellar jitters may prove to 
be a signiflcant source of noise that would rerider it d ifficult to measure sub meter-per-second 
velocity-shifts. More recently, iBouchy et al.l (j2005al ) have actually achieved sub meter-per- 
second accuracy with the HARPS instrument (spectral resolution of 115,000), and they 
have found a dispersion in night-averaged radial velocity measurements for a particular star 
(HD160691) of ~ 0.4 cm s~^ for nights when they took many ( ^ 200) observations. Since 
in our situation (taking spectra during ingress, say) we have minutes, not hours, available, 
the rms scatter in ingress-averaged radial velocity measurements is likely to be larger than 
what they found. In addition to the difficulties posed by several systematic sources of noise, 
achieving sufficient photon statistics will be difficult for two reasons: for a given throughput 
efficiency r], higher spectral resolution means fewer photons per bin; and r] tends to decrease 
with increasing spectral resolution Rs- 

By the mid-1990s, the timeless quest for high-resolution spectrographs reached a mile- 
stone at the Anglo-Australi an Telescope with the development of UHRF and its resolving 
power of up to 1,000,000 (jDiego et al.l Il995l ). Despite impressive throughput relative to 
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previous endeavors, however, its efficiency was insufficient to obtain the sub decameter-per- 
second Doppler precision on a > 7 star that would be req uired for planet searches. With a 
Rs = 600, 000 spectrograph built at Steward Observatory, iGe et al.l (120021 ) obtained stellar 
spectra with Rs ~ 250, 000 and throughput of 0.8%. Furthermore, they predicted t hat by op- 
timizi ng their technology they could increase the throughput to 4%. More recently, iGe et al. 



(120061 ) detected a new planet, around HD 102195, with the Exoplanet Tracker instrument 
at Kitt Peak. This instrument has resolution of Rs ~ 60, 000 and total throughput of 18%. 
Plans for a spectrogra ph that has resolving power of 120,000 on a thirty meter telescope 
(ITokunaga et al.ll2006l ) give cause for optimism that increased aperture area and efficiency 
feeding high and ultrahigh-resolution spectrographs will, in coming years, provide accurate 
enough spectra that tidally locked rotation of HD209458b has a detectable influence. 



3. A Model of a Planetary Transit 

We consider the spectrum of a star whose companion planet transits across the face 
of the stellar disk from Earth's perspective. The primary effect of the planet is to reduce 
the stellar flux at all wavelengths, but the planet's chemical composition, internal structure, 
and rotation rate influence the spectrum in wavelength-dependent ways. Since each of these 
factors - and others too, such as the star's rotation - influences the observed spectrum, we 
built a model that incorporates the many parameters related to each process. The star and 
the planet are both assumed to rotate as solid bodies, with no other (nonthermal) motion in 
their atmospheres. Since deviations from pure solid body rotation are likely to be no more 
than 25% over the disk of the star - e.g . , the Sun's equator-to-pole variation in rotation 



rate is about 21%, as per [Howard et al.l (119841 ). this is probably a reasonable assumption 



for the star. For the planet, this assumption might fail, because wind-speeds in excess of 
the equatorial rotation speed of feq ~ 2 km s~^ are predicted by many models, as described 
in § [T] above. Still, when making this initial study of the spectral effect of the motion of a 
transiting planet's atmosphere, separating rotation from other processes makes the problem 
more tractable. We set parameter values to match measured values from the HD209458b 
system where possible. 

The planet is modeled as an inner component that is entirely opaque and an outer 
component that is isothermal and drops off exponentially. We compute the wavelength- 
dependent optical depth due to the sodium D-doublet at ~ 590 nm in the planet's atmo- 
sphere; important parameters include the temperature and density of the planet's atmosphere 



and its Na-content. We use the Voigt profile - as described by, e.g., iPress fc Rybickil (119931 ) 



- to calculate r[A], the optical depth to absorption along the line of sight. 

As the planet transits the star, there are four points of "contact" between the planet 
and the star (really between their projections on the sky): when the disk of the planet ffist 
touches the disk of the star; when the planet is first entirely over the stellar disk; when the 
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planet is last entirely over the stellar disk; and when the planet last touches the stellar disk. 
We will additionally sometimes refer to "1.5*'^" contact (half-way between first and second 
contact), and analogously to "2.5*'^" and "3.5*'^" contact. 

As described in §|2]above, the type of distortion that a planet's rotation imposes relative 
to a non-rotating planet changes depending on when during the transit the observation is 
made. During ingress or egress, the rotation of a tidally locked planet's atmosphere will 
impose a distortion similar to an overall shift relative to a non-rotating planet: redshift 
during ingress; blueshift during egress. When the planet is in mid-transit, in the middle 
of the stellar disk, the overall distortion to the spectrum imposed by its rotation is akin 
to a star's rotational broadening. Since the line-centers of the lines we are considering are 
heavily saturated and therefore fiat at their cores, rotational broadening has the somewhat 
counterintuitive effect of steepening the cores of the profiles while broadening the wings. 
We will discuss this in greater detail in the next section. Although the tj^e of distortion is 
different during ingress and egress from during mid-transit, it turns out that the amount of 
distortion, in terms of S/N ratio, is nearly constant throughout transit. This, too, we will 
discuss in § m below. 

We simulate the HD209458b system, with a 1.32Rj planet in a 3.5 day orbit, orbiting 
a GO star at with radius 1. 5R(?) th at is 47 pc away. Our model star has the limb darkening 
profile that iKnutson et al. I J2007bh measured for HD209458. In order to approximate the 



fits to the data in ICharbonneau et al.l (120021 ). we assign our model planet's atmosphere a 



sodium-content and cloud deck height (1% solar, and 0.01 bars) that are comparable to 
the parameter-combinations that result in the best fits in that paper. Finally, we present 
results at our simulation's spectral resolution {Rs = 700, 000), and we simulate transit events 
observed using two different lower resolution spectrographs, one with spectral resolution 
R'g = 50, 000 and one with R'g = 150, 000. All spectrographs (and associated optical paths) 
in our simulations have 100% throughput efficiency. In the remainder of this section, we 
provide a detailed description of our parameterization of the problem. 



3.1. Parameters of the Star 

The parameters related to the star are listed in Table [TJ They are set to match mea- 



sured p arameters for IID209458, and we use the limb-darkening profile from iKnutson et al. 



(j2007bl ). We normalize the flux to that of a blackbody of temperature T^, of the size and at 



the distance of HD209458. 
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3.2. Parameters of the Planet 



The parameters related to the planet are in Table [2l We model the planet as an inner 
component that is essentially a billiard ball (completely opaque at all wavelengths) and an 
outer component that is an isothermal atmosphere with scale height H = Rg^sTp/ fig, where 
i?gas is the gas constant, fi is the molar mass, and g is the acceleration of gravity. The 
density of our model planet's atmosphere varies as p = poexp[{r — Rp^)/H], where -Rp^ 
is the radius of the optically thi ck part (some authors have called this radius the "cloud- 
deck" (jCharbonneau et al.ll2002l )). This hypothetical cloud deck could cause the planet to 
be o ptically thick at a highe r alt itude than would othe rwise be expected, as discussed in, 
e.g.. iRichardson et al.l (120031 ) and ISudarsky et al.l (120001 ). The cloud deck causes the optical 
depth as a function of radius in our model to have a singular discontinuity at radius Rp^. 



3.3. Spectral Parameters 

The parameters pertaining to the shape of the observed spectrum are in Table [31 In 
addition to the location of the planet within the stellar disk, the shape of the stellar spectrum 
and the wavelength-dependent opacity of the planet's atmosphere together influence the 
transmission spectrum. 

Spec_Shape is a parameter that can take on the values "Flat", "Blackbody", or "Solar", 
and determines the rest-frame spectrum of the model stellar photosphere. (The integrated 
stellar spectrum is the convolution of the rest-frame spectrum with the stellar rotation 
proflle.) When "Flat" is chosen, the rest-frame model stellar spectrum intensity is set to the 
mean value of the blackbody intensity in the specifled wavelength range [Amin, Amax], which, 
in our simulation, is set to [580 nm, 600 nm] . When "Solar" is chosen, the model stellar 
spectrum instensity is set to a high-resolution solar spectrum that is normalized to the flux 

Table 1. Model Transit Parameters: Star 



Parameter 


Description 


Value 


M, 


Star Mass 


1.05 Mq = 2.09 X 10^3 g 




Star Radius 


1.05 Rq = 7.35 X 10^° cm 


T 


Star Temperature 


6100 K 




Distance to star 


47 pc 




Stellar Rotation Period 


1 month 
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from HD20945qj; but the Na D lines in this high-resolution spectrum have been replaced 
by Gaussian fits to the solar lines. 

The planet's atmosphere has A^^abs absorption features, each of which is due to an element 
with a given fraction of the solar abundance. In the models presented in this paper, N^hs = 2: 
we consider the Na doublet at 588.9950 nm and 589.5924 nm, with sodium at fractional 
abundance = X^a p/^Na. © = 0.01 of the solar abundance. Each line is modeled as a 
Voigt profile, as described in, e.g.. Press & Rybicki (1993). 

3.4. Parameters of Observing and Computing 

The final set of parameters, listed in Table HI includes those that specify the observer 
and those that determine how the observation is discretized for the purpose of numerical 
computation. The model observational setup is determined by three parameters: the tele- 
scope's diameter D (6 m in our simulations) and efficiency t] (100%), and the spectrograph's 
spectral resolution Rs (we set Rs to 700,000 for the purpose of computing the model, and 
we re-bin to lower, more easily achieved resolutions - 150,000 and 50,000 - after computing 
a model). These three parameters prescribe the sizes of the spectral bins and the rate at 
which those bins are capable of collecting light. 

In order to compute the flux at Earth as a function of wavelength, we begin by dividing 
the stellar disk into nf,* concentric annuli, and we divide each annulus into n^^ azimuthal 
sections. In each section, the redshifted spectrum and the normalization must both be 
computed. 

Knowing the stellar rotation rate and axis, we may calculate the recessional velocity 
of any point on the star's surface as a function of its location on the projected stellar disk, 
and we redshift the spectrum from each part of the star accordingly. When the planet is 
in-transit, we separate the stellar disk into an annulus that contains the planet and the rest 
of the disk that we treat as described above. The annulus that contains the planet is treated 
almost as above - divided into Ub^j^ sub-annuli, each of which has ^0^^ azimuthal sections 
- but the sub-annuli are incomplete, interrupted by the planet. 

In order to sample the planet's atmosphere, we divide the region that overlaps the 
star into Ubp concentric annuli around the planet's center, each of which is divided into n^^ 
azimuthal sections. In each section, we must determine the optical depth and multiply by 
exp (— r). In calculating the optical depth, we note that in the case that the planet's rotation 
axis is entirely normal to the line-of-sight, if the planet rotates as a solid body then the 



"'From \ftp://solarch. tuc. noao. edu/\ 
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Table 2. Model Transit Parameters: 


Planet 




A. Cll CUlldjd 




V <X±LIC 




f 


Planet Mass 


0.69Afj : 


= 1.31 X 10^^ g 




Optically Thick Planet Radius 


l.2,2Rj = 


: 9.44 X 10^ km 


Po 


Planet Pressure at Rp^ 




0.01 bars 


H 


Planet Atmosphere Scale Height 




500 km 


T 


Planet Atmosphere Temperature 




1300 K 


/tl 


Prac. Tidal Locked Rot. Rate 


or 1 (Veq = 


or 2 km s"^) 


a 


Semi-Major Axis 




0.046 AU 


#11 


Numl")er of Scale Heights in Atm. 




15 



^Parameter values are set to match measured values from the HD209458b system where 
possible. 



Table 3. Model Transit Parameters: Spectral Features 



Parameter Description Value 



Spec_Shape Shape of Star Spectrum Flat, Blackbody, or Solar 

Amin Min. Wavelength in Sim. 580 nm 

Amax Min. Wavelength in Sim. 600 nm 

-^abs # Abs. Features in P. Atm. 2 

/qi Frac. Solar Abund., First Line 0.01 

Aoi First Line-Center 588.9950 nm 

Akii Transition Prob. First Line 6.16 x 10~^ s~^ 

Qii Stat. Wt. Lower Level First Line 2 

Qki Stat. Wt. Upper Level First Line 4 

/o2 Frac. Solar Abund., Second Line 0.01 

Ao2 Second Line-Center 589.5924 nm 

Aki2 Transition Prob. Second Line 6.14 x 10~^ s~^ 

gi2 Stat. Wt. Lower Level Second Line 2 

gk2 Stat. Wt. Upper Level Second Line 2 



'^In parameters that have i and k subscripts, i indicates the lower level (3si/2 for both 
lines) and k indicates the upper level (3p3/2 for the bluer line and 3pi/2 for the redder line). 
The fractional solar abundance is set to 0.01 in order to achieve modest agreement with data 
observed for the Na D doublet in HD209458b's atmosphere. 
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radial component of its recessional velocity is constant along a ray: 



T[b„<P,\]=N[bp] 



X a 



A 



1 



{vp[bp,(f)p]/c)_ 



where the column density is calculated in terms of a function G that is specified below: 
N[bp] = noG[bp, RpQ, H]. In equation ([8]), Vp[bp,(j)p] is the recessional velocity of the planet, 
as a function of radius and azimuth, which depends upon the orbit and the rotation. Note 
that there is a single Vp along a given hne-of-sight defined by a {bp, (pp) pair only under the 
assumption of solid body rotation. The rest-frame cross-section a[X\ is computed according 
to the Voigt profile. The function G is defined as the following integral: 



G\bp, RpQ, H] 



/oo 



H 



dl 



OO 



bp > Rp 
bp < Rp 



(9) 



4. Model Transit Spectra 



As described in § \2.2\ we seek the expected signal-to-noise ratio for distinguishing 
between the spectrum that would be observed due to a non-rotating planet (or one that is is 
in a Cassini state with its rotation axis nearly in the plane of orbit) and the spectrum that 
would be observed due to a tidally locked planet. The computed model spectrum J\f[\] is 
the time integral of the instantaneous spectrum Af[X\ and consists of the number of photons 
detected per wavelength bin: 

Ar[A]^Ar[A]AUs, 

for some small exposure time Atobs- 

The model signal (of rotation) per bin that we are looking for is the difference between 
the rotating model spectrum A/'rot and the non-rotating model spectrum A/'no rot^ 

Sb = [Kot[X\ - Ao rot[A]) AUs. (10) 

We make the optimistic approximation that the noise per bin is just the photon-noise: 



iVfe= V^norot[A]Atobs. (H) 

The total signal-to-noise ratio in a single exposure, then, is the sum in quadrature of Sb/Nb 
for all wavelength bins A^: 



S/N 



I xv^. (12) 



\ *=1 \ Y-MiorotfA 
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A similar summation in quadrature applies over all exposures. Note that, in principle, 
the expression in equation (IT^ is insensitive to the sizes of bins and hence to the spectral 
resolution Rs, as long as the bins are small relative to the Gaussian width of the absorption 
feature under consideration. Our simulations indicate that the spectral resolution must be 
^ 500, 000 in order for S/N to be nearly independent of Rs- 

The effect of rotation, both during ingress and during mid-transit, is illustrated in 
Figure O For illustrative purposes, in this figure we assume a uniform star (flat spectrum, 
non-rotating, no limb-darkening). On the left panels of Figure [21 we show a snap-shot 
during ingress (at 1.5**^ contact) and on the right, we show a snap-shot during the middle 
of a transit (2.5**^ contact). The quantity plotted is 7?.' = 7?. — 1 from BOl, where 

-^in transit [A, t] ^^^^ 
-^out of transit [-^) ^] 

The bottom panels of Figure [2] show the difference spectra between the models with a tidally 
locked planet and the models with a non-rotating planet (JZrot — T^no rot)- 

As described in §[2] above, a planet's rotation causes the centroids of stellar absorption 
features to shift relative to a non-rotating planet. In Figure [Sj centroid shifts (in velocity 
units) are plotted as a function of position in transit, for a planet transiting in front of a 
realistic star model with a Sun-like spectrum. The recessional velocity increases roughly 
sinusoidally during ingress, reaching a peak of about 60 cm s~^ at l.S*'^ contact. During 
mid-transit, between 2"^*^ and 3'''^ contacts, the net velocity shift is much smaller. Egress is 
nearly perfectly symmetrical with ingress, though the velocity shifts have the opposite sign. 



7^ A 



The cumulative and incremental signal-to-noise ratio across the transit are shown in the 
top and bottom panels, respectively, of Figure HI As a planet proceeds in its transit from first 
contact to fourth contact, the S/N builds up steadily throughout the transit, as can be seen in 
the top panel of this figure. This cumulative increase reflects the steady incremental S/N per 
second of observation, which, for a planet crossing in front of a uniform star, is shown in the 
bottom panel. The three curves in the bottom panel represent, in decreasing order of S/N, 
the incremental S/N curves that are expected for a spectrograph with ultrahigh-resolution 
(the simulation's resolution Rs = 700,000), for a spectrograph with as high-resolution as 
bHROS on Gemini [R'g = 150,000), and for a spectrograph with more standard high- 
resolution {R'g = 50,000). In the top panel, we show the cumulative S/N (assuming the 
simulation's resolution) for a planet in front of a uniform star, and for a planet in front of a 
more realistic sun-like star that rotates once per month, has a limb-darkening profile, and 
has a solar spectrum. It is apparent in the top panel of Figure H] that a spectrograph with our 
simulation's resolution would see the effect of rotation at S/N = 7.1?7^/^ in a single transit 
in the case of the simplified uniform star, and at S/N = 3.2r]^^^ in one transit of the realistic 
star. The effect of including in the simulation the realistic features of stellar rotation, limb- 
darkening, and a solar spectrum is therefore to depress the S/N of the effect of tidally locked 
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Fig. 2. — Upper panels show snap-shot spectra for one of the the Na D hnes for two different 
model planets (tidally locked and non-rotating); lower panels show the difference between 
the two model spectra. The quantities plotted are TZ' = 71 — 1 (upper panels) and ATI' 
(lower panels), where 7l[X,t\ = iVin transit [A, t]/iVout of transit [A, t]. In the upper panels, the 
blue curve is the tidally locked planet's transit spectrum, and the green curve is the non- 
rotating planet's transit spectrum. In the lower panels, the difference between the rotating 
and non-rotating planet's spectra. Left: Halfway through ingress (at l.S**^ contact). Right 
Halfway through the whole transit (2.5*^ contact). 

rotation by a factor of slightly more than 2. The bottom panel of the figure indicates that 
our predicted S/N for one transit of a realistic star (3. 2?]^/^) might be adjusted downward 
by ~ 50%, depending on the spectral resolution, indicating a total one-transit S/N for the 
case of the realistic star of about ~ l.Tr/^/^. Finally, we do not show the incremental S/N for 
the realistic star model, but, as the top panel of the figure indicates, a larger proportion of 
the total signal comes from ingress and egress in the realistic star model than in the uniform 
star model. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Fraction from First Contact to Fourth Contact 



Fig. 3. — Centroid-shift of Na D lines from tidally locked rotation from the beginning 
to the end of a transit, relative to an identical but non-rotating planet; Sun-like stellar 
spectrum. The vertical lines denote l.S*'^ and 3.5*^ contact (cyan) and second and third 
contact (magenta). Between first and second contact, the spectrum with the rotating planet 
is redshifted relative to the non-rotating planet by up to about 60 cm s~^; between third 
and fourth contact, it is blueshifted by the same amount. This plot samples the transit at 
60 regularly-spaced points. Parameters were chosen to represent the HD209458 system. 

5. Discussion 

The S/N scales with the square root of the number of absorption lines under consider- 
ation, with the square root of the number of transits, and with several other parameters, as 
follows: 




i^P[R's])iS/N), (14) 

In this equation, Ai'abs is the number of absorption lines, D is the telescope's diameter, t] is the 
efficiency of the instrumental setup, K is the star's magnitude, Vq is HD209458's magnitude 
(about 7.6), is the stellar radius, Rq is HD209458's radius, Rp and H are the planet's 
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Fig. 4. — Top: S/N for distinguishing between a rotating and a non-rotating planet, accu- 
mulated during a single transit. Upper curve displays the cumulative S/N for the distinction 
of a tidally locked planet from a non-rotating model, assuming a completely uniform star: 
the model star has a flat spectrum, is not rotating, and there is no limb-darkening. Lower 
curve displays the same, but for a more realisti c star: the model star r otates once per month, 
has the limb-darkening profile determined by lKnutson et al.l (j2007bl ). and has a solar spec- 
trum. These curves assume the spectral resolution used in the simulation {Rs = 700,000); 
for more easily achievable spectral resolution, these curves are adjusted downward by 30% 
- 80%, as displayed in the bottom panel. The vertical lines denote 1.5**^ and 3.5*^^ contact 
(cyan) and second and third contact (magenta). Also, both top and bottom panels show 
model results assuming 100% throughput efliciency (r/ = 100%). In practice, these curves 
would be adjusted downward by a factor i]^^"^, or ~ 

0.04^/2 ^ 0.2 for an ultrahigh reso- 
lution spectrograph. Bottom: Incremental S/N for sixty 1-second observations across the 
transit, assuming a uniform, non-rotating star with a flat spectrum. Top curve shows the 
incremental S/N, assuming a spectrograph with the spectral resolution of the simulation 
(i?5 = 700,000). Middle and bottom curves show the same, assuming spectral resolutions 
R'g = 150, 000 and R'g = 50, 000, respectively. Parameters of the system (except for the 
stellar spectrum in the bottom panel) were chosen to represent the HD209458 system. 
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radius and scale height, and and Hq are HD209458b's radius and scale height. 'iplR's] 
is a function that characterizes the (nonlinear) response of S/N to the spectral resolution 
of the spectrograph; for instance, ^^[700, 000] ^ 1, ^[150,000] ~ 0.6 and ^[50,000] ~ 0.2. 
Finally, {S/N)q = 7.1 for a uniform star and p/^TSSa © = 0.01; and {S/N)o = 3.2 for a 
more realistic (rotating, limb-darkened) star with a solar spectrum. 

Culling the signal from mid-transit, however, is much more difficult than from ingress 
and egress, because the shape of the distortion depends more sensitively on the structure of 
the planet. For the realistic star, our model predicts S/N ~ 1.3?]^/^ (with the 6 m telescope, 
with Poisson-dominated photon noise) for ingress alone, and therefore a factor 2^/^ greater, 
or S/N l.8r]^^^ for ingress and egress. 

In Table El we present the number of orbits that must be observed in order to make a 
5a detection of rotation, for various combinations of parameters. In all cases, we assume a 
realistically achievable optical setup, with spectral resolution R'g = 150, 000 and throughput 
1] = 4%. Improvement in R'g will yield a modest improvement in S/N, and improvement 
in rj could be quite significant. In the third column of Table [5] (A), we present the total 
available S/N for the whole transit, according to our model, while in the fourth column (B), 
we present the available S/N from just ingress and egress. The fifth and sixth columns are 
based on the S/N in the fourth column. 

Finally, as a sanity ch eck, our model can be tested by comparing it to the analysis of 



Charbonneau et al.l (120021 ). That analysis suggests that the sodium content of a planet's 
atmosphere can be determined by comparing the flux in a narrow band centered on the 
sodium resonance lines with the flux in a wider band surrounding but excluding the lines. 
In that paper, the decrement in the narrow band containing the sodium features is named 
nNa, and was measured to be —2.32 x 10"'' for HD209458 during the middle of transit. They 
presented several models, all of which over-predicted the sodium decrement. The model 
that predicted the smallest magnitude of the decrement had 1% solar metalicity and cloud 
deck at 0.0368 bar; this model predicts nNa ~ —3.4 x 10~^ in mid-transit. Our model (1% 
solar metalicity, cloud deck at 0.01 bar) predicts ~ —4.1 x 10~^ in mid-transit, as shown in 
Figure [51 We conclude that our model is in reasonable agreement with both Charbonneau 
et al.'s model (~ 20%) and with the actual data (~ 40%). 



6. Conclusion 

Our investigation indicates that, with currently available instruments, it will be difficult 
to obtain the sensitivity needed to achieve a minimal S/N 5 detection of tidally locked 
rotation of the planet HD209458b. Nevertheless, it appears that the effect of rotation will 
have significant (though small) infiuence on transit-spectra taken with current or near-future 
instruments. Because this infiuence is so small, it is worth considering ways that the S/N 
could be improved, according to the scalings summarized in equation (IT^ . 



Table 4. Model Transit Parameters: Observation and Computing 



Parameter Description Value 



D Telescope Diameter 2.4 m — 30 m 

T] Spectroscope Efficiency 1.00 

R's Obs. Spec. Resolution 50,000 - 700,000 

Tint Integration Time 932.088 s 



Rs Comp. Spec. Resolution 700,000 

At Time-Step in Integration 50 s 

rib* # of Star Annuli 10 

n^^ # of Star Azimuthal Sections 16 

nb^:A 7^ of S. Annuli in P. Annulus 10 

n^^^ # of S. Azim. See's in P. Ann. 10 

Ubp of Planet Atm. Annuli 20 

n^p # of Planet Atm. Azim. Sections 20 



*^Parameter values are set to match measured values from the IID209458b system where 
possible. 



Table 5. Required Number of Transits For 5cr Detection 



Spec_Shape 


D 


(S/N)/Tr. (A) 


(S/N)/Tr. (B) 


Req. # Transits (A) 


Duration (A) 


Flat (uniform star) 


6 m 


1.1 


0.52 


~ 25 


~ 3 months 


Solar 


6 m 


0.48 


0.27 


~ 110 


~ 1.25 years 


Solar 


10 m 


0.79 


0.45 


~ 40 


~ 5 months 


Solar 


30 m 


2.4 


1.3 


~ 5 


~ 3 weeks 



^A table showing the number of orbits required to achieve a. S/N = 5 detection of tidally 
locked rotation of HD209458b, assuming a planetary Na-content that is 1% solar. As- 
sumptions: we are monitoring two spectral lines (the Na doublet) with a spectrograph with 
resolving power R'g = 150, 000, and an efficiency in the optical setup of 0.04. The third col- 
umn "(S/N)/Tr. (A)" is the S/N per transit, as shown in Figure IH and the fourth column 
is the S/N per transit, considering only ingress and egress. The "(A)" in the fifth and sixth 
columns indicates that the numbers of transits and the required duration are tabulated for 
S/N values from column three. 
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Fig. 5. — The flux decr ement due to sodium in th e planet's atmosphere, as shown by the 
quantity n^^ defined by ICharbonneau et al.l (120021 ) . This decrement reaches a maximum in 
magnitude of —4.1 x 10~^, halfway through transit, in our model. 



The most obvious way to boost S/N is to increase the collecting area of the telescope, 
the effect of which is shown in Table |5l With an extremely large (30 m-class) tele- 
scope, the required discrimination power could be achieved in ~ 5 — 15 transits, if our 
optimistic assumptions about the noise are not too far off. 

As the abundances of other elements or molecules in the planet's atmosphere are iden- 
tified, the number of absorption lines A^abs that can be utilized can increase. Observing 
four absorption features instead of two will boost the S/N per orbit by 2^^^ and will 
therefore cut in half the required observing time to achieve a fixed target S/N. 



The HD189733 system is only 19.3 pc away (iBouchy et al.l l2005bl ). but the star is 
smaller and dimmer than HD209458, and has the same apparent magnitude. We find 
that, overall, it has no relative advantage in terms of S/N. It may be unlikely that 
a star significantly brighter than HD209458 will be found to have a transiting giant 
planet companion, but if this should happen the signature of rotation will be more 
readily apparent in the system with the brighter star. 

Finally, since S/N is proportional to the size of the atmosphere, represented by S^tm 
in equation ([7]), a planet with a more extended atmosphere would show the distortion 
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due to rotation more readily. The size of the atmosphere depends on the product of Rp 
and the height of the atmosphere Ilatm ~ 15H, so the S/N could be increased in the 
case of a larger planet, a hotter one, or one with lower surface gravity. Unfortunately, 
although some other planets are somewhat larger, none of the puffier planets that have 
been discovered so far happen to transit a star as bright as HD209458. If a more 
bloated planet is found to transit a star as bright as HD209458, such a planet would 
probably be a more promising target for the type of study described in this paper. 



Although the optical transit spectrum of an extrasolar giant planet contains many pre- 
cious clues regarding the nature of rotation and climate or weather, we emphasize that we 
have so far addressed only the forward problem: deep, high spectral resolution transit spec- 
tra will be distorted by HD209458b's tidally locked rotation, if indeed it is tidally locked. 
But seeing the effect of rotation in a spectrum and knowing that we are seeing the effect of 
rotation are two very different things. The inverse problem - identifying the rotation rate 

- will be much more challenging, and is beyond the scope of this paper's investigation. In 
brief, it is impossible to decode the clues that the spectrum holds without a Rosetta stone 

- namely, without already possessing accurate knowledge of abundances in the planet's at- 
mosphere. We note that when it comes to solving the inverse problem, it may be possible to 
make do without a perfect Rosetta stone by comparing the star's spectrum at two different 
times of transit instead of comparing it to an accurate input model. Specifically, comparing 
the ingress spectrum to the egress spectrum should eliminate various systematic uncertain- 
ties that may effect both ingress and egress spectra. This strategy involves foregoing spectra 
between second and third contacts, at a cost of a factor of ~ 2 in S/N, but has the ad- 
vantage of being more robust with respect to error s in modeling the planetar y atmosphere. 



Furthermore, if detailed analysis along the lines of ICharbonneau et al.l (|2002| ) should reveal 
many more planetary absorption features, a large current-generation telescope could detect 
the motion of an extrasolar planet's atmosphere in of order a month. 

A final point to keep in mind is that the technique presented in this paper is not limited 
to discerning the difference between a tidally locked planet and a non-rotating planet; this 
technique can be used to investigate the applicability of any model of the motion of a 
transiting planet's atmosphere. Some atmospheric circulation models of Hot Jupiters predict 
extremely fast winds, up to ~ 9 km s~^, or more then four times the predicted equatorial 
rotation rate. It might therefore turn out to be easier to detect winds than to detect rotation. 
Moreover, if an analysis similar to the one described here were to find a deviation from the 
hypothesis of tidally locked rotation, this would mean that something interesting is going on 
- the planet could be rotating roughly as a sohd body but at a different rate; there could be 
significant non-rotational motion in the atmosphere; or there could be some combination of 
these phenomena. 
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